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ABSTRACT 



I Context. The inverse Compton catastrophe is the dramatic rise in the luminosity of inverse-Compton scattered photons predicted 
^ to occur when the synchrotron brightness temperature exceeds a threshold value, usually estimated to be lO^'^ K. However, this 
efltect appears to be in contradiction with observation because: (i) the threshold is substantially exceeded by several intra-day 
_ variable radio sources, but the inverse Compton emission is not observed, (ii) powerful, extra-galactic radio sources of known 
\^ ■ angular size do not appear to congregate close to the predicted maximum brightness temperature. 

T-H ' Aims. We re-examine the parameter space available to synchrotron sources using a non-standard electron distribution, in order 
to see whether the revised threshold temperature is consistent with the data. 
' Methods. We apply the theory of synchrotron radiation to a population of monoenergetic electrons. The electron distribution 
' and the population of each generation of scattered photons are computed using spatially averaged equations. The results are 
, formulated in terms of the electron Lorentz factors that characterise sources at the threshold temperature and sources in which 
^ I the particle and magnetic field energy density are in equipartition. 

' Results. We confirm our previous finding that intrinsic brightness temperatures Tb ^ 10^'* K can occur without catastrophic 
, ^ cooling. We show that substantially higher temperatures cannot be achieved either in transitory solutions or in solutions that 
' balance losses with a powerful acceleration mechanism. Depending on the observing frequency, we find strong cooling can set 
, in at a range of threshold temperatures and the imposition of the additional constraint of equipartition between particle and 
• • ■ magnetic field energy is not warranted by the data. 
. ■ Conclusions. Postulating a monoenergetic electron distribution, which approximates one that is truncated below a certain Lorentz 
, factor (7min), alleviates several theoretical difficulties associated with the inverse Compton catastrophe, including anomalously 
5_j high brightness temperatures and the apparent lack of clustering of powerful sources at lO'^^ K. 
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1. Introduction 

When interpreting the synchrotron spectrum of power- 
ful extra-galactic radio sources, it is usual to assume 
that the underlying electron distribution has a power- 
law form dA^/d7 oc 7""^. For a homogeneous source, the 
spectrum peaks at the point v — z^abs, where the opti- 
cal depth to synchrotron self-absorption is of the order 
of unity, and above this frequency the intensity falls off 
as Ii, cx i^^f''"^)/-^. Because more electrons become ef- 
fective at absorbing the radiation as the frequency de- 
creases, the optically thick part of the spectrum is not of 
the Rayleigh- Jeans type, but has instead I^, cx v^^^, inde- 
pendent of the power-law index of the underlying distri- 
bution, (provided q > 1/3). Correspondingly, the bright- 
ness temperature peaks at « t'abs, falling off as z^^/^ to 
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lower and as ly ('?+^)/^ to higher fr equencies. In this case, 
Kellermann fc Paulinv-Tothi (|l969f ) found that the ratio 
of the luminosity Lie carried away by inverse Compton 
scattered photons to that of synchrotron photons is 



Lie 



thresh 



Tr 



thresh 



(1) 



where Tr is the intrinsic brightness temperature at the 
peak of the radio emission and Tthresh ~ 10^^ K, depending 
somewhat on the parameter q and the maximum frequency 
at which synchrotron radiation is emitted, corresponding 
to an assumed cu t-off in the power-law electron spectrum 
()Readheadlll99l . 

The rapid rise in total luminosity implied by Eq. ([1]) 
when Tr exceeds Tthresh is called the "inverse Compton 
catastrophe" . The implied energy requirement of a IJy 
source at z = 1 is ~ lO'^^ergss"^ at Tr = lO^^K, 
and a prohibitive ~ lO^^ergss"^ at Tr — 10^^ K at 



2 



O. Tsang and J. G. Kirk: The inverse Compton catastrophe 



IGHz. A source of 1 GHz photons with magnetic field 
~ 1 fiG, boosts a radio photon into the X-ray band in a 
single inverse Compton scattering. At Tb = 10^^ K this 
would imply a nanoJansky X-ray flux, not untypical of 
strong extragalactic sources. However, at Tb = 10^^ K 
Eq. H]) predicts an X-ray flux at the milliJansky 
level, in contradiction with o bservation (Fossati et a l.. 

200d: iTayecchio et all 



1998t iSambruna et al.l 

Padovani et al.l 12004 ICuainazzi et al.ll2006[ ). 



Despite this, several sources that display intra-day 
variability (IDV) in their radio emission have an im- 
plied brightness temperature that exceeds 10^^ K by sev- 
eral orders of magnitude, i f the observed variability is 



intrinsic (e.g., iKraus et~al 
ity is caused 



20031 ). Even if the variabil- 
by scintillation, the implied brightness 
tempera ture can still greatly exceed 10^^ K for some 
sources (IWagner fc Witzellll995h . Currently, the most ex- 
treme example is the source PKS 0405-385. This source 
displ ays diffractive scintillation (jMacquart &: de Bruvn 
20051) . which places an upper limit on its angular size that 
corresponds to a brightness temperature of 2 x 10^"* K. 
These sources are generally assumed to be beamed, i. e., to 
be in relativistic motion t oward s the observer (e.g., Reea 
1966[ IJones fc Burbidgd Il973t ISingal k Gopal -Krishna' 
19851 ) . In this case the intrinsic temperature is lower than 
that deduced for a stationary source by a factor of T)'^ (for 
intrinsic variability) or V (for scintillating sources). (Here 
the Doppler factor T> = y^l^^fp / (1 — /3 cos 9) with /3c the 
source velocity and 9 the angle between this velocity and 
the line of sight.) Nevertheless, the observed brightness 
temperatures are too high to be accounted for by Doppler 
factors similar to th ose estimated froni observations of su- 
per luminal motion ( Cohen et al. 20031) . 



A second problem arises with powerful sources whose 
angular extent can be measured directly. In an analysis 
of high brightness temperature radio sources in which 
Doppl er beaming is thought to be absent, iReadhead 
(|l994[ ) measured a brightness distribution that cuts off 
at 10^^ K; one order of magnitude lower than the in- 
verse Compton limit. This appears consistent with ob- 
servations of a sample of 48 so urces showing super lumi- 



nal motion (jCohen et al.l 120031 ) , in which it was found 



that the intr insic bri ghtness temperature cluster around 
2 X 10^° K. IReadhead (1994) argued that an apparent 
maximum brightness temperature significantly lower than 
10^^ K could not be caused by catastrophic Compton cool- 
ing. Instead, he suggested that sources arc driven towards 
equipartition between their magnetic and particle energy 
contents. Assuming, in addition, that observations are 
taken at the peak of the synchrotron spectrum, and that 
the electron distribution is a power-law, he showed that 
the expected distribution of brightness temperatures was 
consistent with that observed. 

In this paper, we re-examine these two problems as- 
suming that the source contains a monoenergetic elec- 
tron distribution instead of the conventional power-law. 
Although this assumption appears at first sight highly re- 



strictive, the form of the synchrotron emissivity means 
that under some circumstances such a distribution pro- 
vides a good approximation to several more commonly 
encountered cases, including that of a conventional power- 
law distribution that is truncated to lower energy at a 
Lorentz factor 7min- Such distributions have been pro- 
posed in connection with radio sources for a variety of 
reasons: the absence of low energy electrons can account 
for the lack of Faraday depolarisation in parsec- scale 
emission regions (jWardlc 19771 : IJones fc Odel]||l977n and 
has recently been discussed in connection with statistical 
trends in the observed distribution of superluminal ve- 
l ocities as a function of obs erving frequency and redshift 
(jGopal-Krishna et al. I l2004l ). Also, iBlundell et all (|2006l ) 
recently examined the radio and x-ray emission from the 
lobe regions of a giant radio galaxies 6C 0905-1-3955, and 
deduced a low energy cutoff of the relativistic particles in 
the hotspots of 7i„in ~ 10^. 

In Sect. [2] we use standard theory to discuss the gen- 
eral properties of the synchrotron spectra emitted by a 
homogeneous source. A set of spatially averaged equa- 
tions describing the evolution of the electron Lorentz 
factor and both the synchrotron and the associated in- 
verse Compton scattered emission is presented in Sect. |3l 
Having identified in these equations the threshold for the 
inverse Compton catastrophe, we discuss the parameter 
space available to stationary s olutions in Sect.lH H ere we 
confirm the results reported in lKirk fc Tsangl ( 20061 ). indi- 
cating that temperatures considerably in excess of 10^^ K 
are permitted. We also show that in the case of resolved 
sources, the onset of catastrophic cooling occurs over a 
wide range of temperatures, consistent with the observed 
temperatur e range . Fina lly, we address in Sect.[5]the sug- 
gestions by ISlvshI (|l992l ) that extremely high brightness 
temperatures can be achieved in nonstationary sources ei- 
ther by injecting electrons at high energy, or by balancing 
their cooling against a powerful acceleration mechanism. 
A summary of our conclusions is presented in Sect.[6l 

2. Synchrotron spectra 

We consider a homogeneous source region characterised 
by a single spatial scale R, that contains monoener- 
getic electrons and possibly positrons of Lorentz factor 
7 and number density rie immersed in a magnetic field 
B. Expressions for the synchrotron emissivity and ab- 
sorpt ion coefficients can be found in many excell ent texts 
(e.g., Rvbicki fc Lightnianl ( 19791 Chapter 6), and Longair 



(|1992L chapter 18)) and are summarised in our notation 
in Appendix lAl 

For any given source there exists a frequency i/abs 
below which absorption is important. Since B and 7 
also define a characteristic synchrotron frequency i^g (see 
Eq.O]), the sources we consider can be divided into two 
categories: those with weak absorption in which z^abs < t's 
and those with strong absorption I'abs > i^s- Note that this 
division is independent of the observing frequency, since it 
relates only to intrinsic source properties. The synchrotron 
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Fig. 1. The synchrotron spectra (upper panel) and 
brightness temperatures (lower panel) of sources with mo- 
noenergetic electrons in the case of strong (blue) and 
weak (red) absorption. The green curves show the opti- 
cally thick [ly = Sv, see Appendix \K\ and optically thin 
{Ijj — TgS^) approximation. In the upper panel, is in 
arbitrary units, and in the lower, the brightness tempera- 
ture is normalised to the energy of the electron. The ab- 
scissa X is the ratio of the frequency to the characteris- 
tic synchrotron frequency of the electrons Vs- The blue 
(red) curves correspond to a source which has an optical 
depth of unity to synchrotron self-absorption at a: w 5 
{x « 0.05). For ease of display, the upper panel compares 
sources with equal flux at high frequency, whereas the 
lower compares sources with equal flux at low frequency. 

spectra that emerge in these two cases are quite different, 
and are illustrated in Fig. [TJ A feature they have in com- 
mon is that the low energy spectrum has the Rayleigh- 
Jeans form cx v"^ , where 1^ is the speciflc intensity at 
frequency v. This property contrasts with the v^/'^ de- 
pendence of at low frequencies of a source containing 
a power-law distribution of electrons. The reason is that 
a power-law distribution contains cold (low energy) elec- 
trons that contribute to the absorption at low frequencies. 

The brightness temperature, Tb = c^IuH'^-v^k-Q), 
where is Boltzmann's constant, is a function of fre- 
quency and is also illustrated in Fig. [TJ At low frequency, 
it attains its maximum value roughly in "equilibrium" 
with the electrons: 7b, max = i^mc^ / ^k^, then decreases 
monotonically to higher frequencies. In the case of weak 



absorption, Te^max oc v^^^^ for I'abs < v < Vs, and then 
cuts off exponentially as z^~^/^exp {—v/vs) once Vg is ex- 
ceeded. In strongly absorbed sources, the brightness tem- 
perature remains almost constant until the frequency ex- 
ceeds Vs upon which it falls off as until the source 
becomes optically thin, after which the exponential cut- 
off T cx v~^/'^ exp {—v/vs) takes over. 

Although four parameters (7, rtc, B and R) are needed 
to define a source model, the division between strong and 
weak absorption is simple. It occurs at a critical Lorentz 
factor 7c given by (see Eq. lA.lip 



7c = 324 x(-^) 
v 1 cm / 

or, equivalently, 



1/5 



R 



7c = 4451 X T^^^ 



B 



ImG 



1 kpc 

-1/5 



1/5 



B 



ImG 



-1/5 



(2) 
(3) 



where tt = UcRcft is the Thomson optical depth of the 
source. Strong absorption occurs for low Lorentz factors 
7 = 7/7C < 1 and weak absorption for high Lorentz fac- 
tors 7 > 1. If the Lorentz factor 7 is held constant, the 
strong absorption regime may be reached from the weak 
by increasing tt at constant i?, or by decreasing B at 
constant tt. 

In his m odel of high-brightness temperature sources. 



SlyshI (jl992f ) considered the strong absorption case. The 
most important property of the assumed distribution in 
this case is the lack of high energy electrons: the addition 
of a population of cold electrons, which would correspond 
to a power-law distribution truncated to higher Lorentz 
factors, would reduce the brightness temperature of the 
source at x < 1 (in Fig. [1]) but would not significantly 
influence this quantity, for x > 1. 

On the other hand, I Crusius- Waetzell (|l99ll ) and 



Prothero^ ( 2003 ) considered weak absorption, where the 
key property of the model distribution is the absence 
of low energy electrons. In this case, the monoenergetic 
model is a good approximation to a power-law distribu- 
tion truncated to lower electron energies at 7 = 7min. The 
addition of a high-energy power-law tail affects the spec- 
trum at a; > 1, but does not change the maximum bright- 
ness temperature achieved at a; <^ 1. Furthermore, the 
truncation need not be sharp: provided the opacity at low 
frequencies is dominated by the contribution of electrons 
with 7 « 7min, the monoenergetic approximation is good. 
This is the case if, for 7 < 7min, the spectrum is suffi- 
ciently hard: diV/d7 cx 7^"? with q < 1/3. In particular, 
the low energy tail of a relativistic Maxwellian distribution 
{q = —2) falls into this category. 

In contrast to the pure power-law distribution, where 
the self-absorption turnover is strongly peaked, the emis- 
sion of a weakly absorbed source — shown in red in the 
upper panel of Fig. [1] — is flat over nearly two decades in 
frequency. It therefore provides a natural explanation of 
compact flat-spectrum sources, eliminating the need to ap- 
peal to a "cosmic conspiracy" behind the superposition of 
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peaked spectra from d ifferent parts of an inhomogeneous 
source ( Marscher Il980l) . 

For the treatment of inverse Compton scattering, it is 
necessary to evaluate the the energy density Us in syn- 
chrotron photons in a given source. To do this, must 
be integrated over angles and over frequency. The result 
depends on the geometry and optical depth as well as the 
position within the source. An average value can be es- 
timated by introducing a geometry dependent factor, C, 
defined according to: 



Us 



AttC 



diy (/.) 



(4) 



where (7^) is conveniently taken to be the specific intensity 
along a ray path that is within the source for a distance R 
and is perpendicular to the local magnetic field. Prothero3 
(l2002h lias evaluated C for several interesting special cases. 
For a roughly spherical source, it is of the order of unity. 
Below, we show that the choice C = 2/3 is consistent with 
our spatially averaged treatment of the kinetic equations. 
The dominant contribution to the integral over the spec- 
trum arises from photons of frequency close to i^s in the 
case of weak absorption, and close to Vabs in the case of 
strong absorption. Using this approximation, for weak ab- 
sorption (7 > 1): 



Us 



4.1 X 10^S\ 



1 cm 



R 



1 kpc 



or, equivalently. 



Us - 2YttC 



VStt 



and for strong absorption (7 < 1): 



Us « 8.9 X 10-^«7jC (f^ 



B 



ImG 



(In 7)^ 



(5) 
(6) 

(7) 



An approximation that is accurate for all values of the 
optical depth is given in Eq. (jA.24p . 

3. Spatially averaged equations 

An approximate, spatially averaged set of equations gov- 
erning the energy balance of particles and synchrotron ra- 
di ation in a source can be foun d fo llowing the approach 
of Lightman fc Zdziarskil (|l987l) and iMastichiadis fc Kirk 
( 1995f ). In terms of the time-dependent synchrotron radi- 



ation energy density Uo{t) one can write: 

^ + c{a,)Uo + ^Uo - (>> 



(8) 



The second and third terms on the left-hand side of this 
equation represent the rate of energy loss by the radia- 
tion field due to synchrotron self-absorption and escape 
through the source boundaries; the right-hand side is the 
power put into radiation by the particles. The angle brack- 
ets indicate a frequency and angle average, but, within this 
spatially- averaged treatment, an exact calculation of the 
frequency average is unnecessary; it suffices to replace the 



absorption coefficient by its value where the energy den- 
sity of the synchrotron spectrum peaks i.e., aX v — in 
the case of weak absorption and v = v^hs in the case of 
strong absorption. In terms of the optical depth to syn- 
chrotron absorption at this point, Tp < 1, the equation 
becomes: 



dC^ 
dt 



c 

'r 



1 



(7o 



0.) 



(9) 



The right-hand side of this expression can now be found 
by demanding it gives the correct steady solution at both 
large and small optical depth. The resulting equation is: 



dt 



R 



(1 



[c/ow - Usm 







(10) 



where Us is the steady-state synchrotron radiation energy 
density, evaluated according to Eq. (jH) , with an appropri- 
ate value of the parameter C,. 

The corresponding equation for the particles that takes 
into account synchrotron absorption and emission as well 
as an acceleration term takes the form 



2^7 

ricrac — — 
dt 



c c , 



(11) 



The first term on the right-hand side of Eq. (|lip is the 
power taken from the radiation field by self-absorption and 
the second term is that returned to it — both of these ap- 
pear in Eq. (jlO[) . The third term describes the energy input 
by particle acceleration. The particular scaling used fol- 
lows that of Slvshl ( 1992h . and models a generic first-order 
Fermi process. For a independent of 7, the acceleration 
rate is proportional to the gyro frequency, and for a = 1 
it equals this value. The acceleration timescale equals the 
crossing time of the source when a = jmc^ /{eBR). 

Multiple inverse Compton scatterings can be ac- 
counted for as follows: First we label the photons present 
in the source according to how many scattering events 
they have suffered after production by the synchrotron 
process. The energy density of these photons is denoted 
by Ui Thus, i — corresponds to photons emitted by 
the synchrotron process which have not undergone a scat- 
tering, and the corresponding energy density is governed 
by Eq. (fTO|) . Assuming the source is optically thin to 
Thomson (or Compton) scattering, the dominant loss 
mechanism for the energy density of photons belonging 
to a given generation i > 1 is escape from the source, 
rather than conversion to the i + I'th generation. In this 
case, we can write for the time-dependence of Uf. 



m 

dt 



R 



:U, 



(12) 



where Qi is the rate per unit volume at which energy is 
transferred into photons of the i'th generation by inverse 
Compton scattering, for i > 1, or by synchrotron radiation 
for i = 0. 

If the inverse scattering process proceeds in the 
Thomson regime a simple expression can be found for Qi. 
However, as i increases, hvi also increases, eventually be- 
coming comparable to the electron energy when viewed in 
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its rest frame. When this happens, Klein-Nishina modifi- 
cations to the Thomson cross section become important, 
reducing the value Qi. We take approximate account of 
this effect by limiting the number of scatterings to iVmax, 
and using the Thomson approximation to evaluate Qi for 
i < Nma,x- In this case, the average energy of a scattered 
photon of the i'th generation is Vi = Aj'^i^i^i/S and the 
rate of such scatterings in unit volume of the source is 
rLcaTcUi-i/ {hvi^i). Therefore 



Q^ 



icUi-i/R for 1 < i < iV„ 



for i > N„ 



(13) 



where the parameter ^ is defined as 

47^ tt 



(14) 



The appropriate value of A^max is chosen by requir- 
ing the average energy of the iVmax generation of photons 
viewed in the electron rest frame 7(47^/3)^°""=/ii'o to be 
less than the electron energy: 



-/Vmax = floor 



In (mc^/ft.t'o) 
21n7 



1 



For synchrotron radiation, Eq. (jlOp implies 

Qo = ^ (Us - Uo) + 



R 



R 



(15) 



(16) 



In the stationary case, Uq — Ug, Eqs. and give 
Qi/Qo = ^- However, assuming scattering in the Thomson 
regime, the ratio of the energy lost by synchrotron radia- 
tion to that by inverse Compton scattering in the steady 
state equals the ratio of the energy density of the magnetic 
field to that of the target photons Qo/Qi ~ i3^/(87rC/i_i), 
which, for i = 0, implies Us = C(i3^/87r). Comparison 
with Eq. ^ then confirms that the spatially averaged ki- 
netic equations are consistent with the choice ^ = 2/3 for 
the geometry dependent factor. Finally, the electron equa- 
tion (jlip acquires the additional loss terms from inverse 
Compton scattering: 



2d7 \ - ^ ^ 

neinc — — = — > Qi + a enc 
at 



(17) 



The set of equations and p?|) can be rewritten by 
introducing the total energy density of scattered radiation: 



N„ 



Ut = ^ C/. 



(18) 



Then, using dimensionless variables according to U ■ 
U {Stt/B"^), i = tc/R and = 8ttcQi/{RB'^) one finds 



dC/x 
"dT 



[1 - t>T = e (■ 



Uo - Un„ 



(19) 



If Um^^^ remains always negligibly small, then all signifi- 
cant scatterings occur in the Thomson regime, and the set 
of equations (fT2|) (for z = 0), (fT7|) . and (fT9|) can be conve- 
niently formulated in terms of three characteristic values 
of the Lorentz factor: 



d£ 
d7 



1 - (7/7cat)^ Ut -t- (7/7cat)^ Uq 



-Uo + Qo 



dF ^ 



Qa + (7/7cat) Ut 



(20) 
(21) 
(22) 



where 7eq is chosen so that there is equipartition between 
particle and magnetic energy densities for 7 = 7oq: 



7cq = / {STTricmc^) 

7cat is given by setting $ = 1 



7cat 



3 

4tt 



(23) 



(24) 



and 7tr corresponds to the maximum Lorentz factor of 
a particle that can be confined in the source, i.e., whose 
gyro-radius is less than R: 



7tr 



eBR/imc^) 



(25) 



The significance of 7cat can be seen from the 
steady state solution of Eqs. (PU)) and (PT|) : Ut = 
Us/ (7cat/7^ ~ 1) • Foi' values of 7 that approach 7cat from 
below, the energy density in the radiation field, and, 
hence, the luminosity diverge. Thus, under the assump- 
tion that all scatterings take place in the Thomson limit, 
no stationary solutions can be found for 



7 > 7cat 



(26) 



This phenomenon is the nonrclativistic or "Thomson" 
manifestation of the Compton catastrophe described in the 
Introduction. In the weak absorption limit. Us = 7^/7catJ 
confirming the well-known result that the Compton catas- 
trophe sets in when the energy density in synchrotron pho- 
tons exceeds the magnetic energy density. However, this 
result does not apply to the case of strong absorption, 
where we find Us ~ 7cat/7c ^ 1- I^ this regime, the syn- 
chrotron radiation energy density can be much smaller 
than the energy density in the magnetic field at the point 
where catastrophic cooling sets in. Physically, the scat- 
tered photons feed on each other to produce the catas- 
trophe in this regime, and do not require a substantial 
synchrotron photon density. In a realistic model, the di- 
vergence of the luminosity is prevented by Klein-Nishina 
effects, that effectively truncate the series in Eq. (IT5|) . For 
example, if Te^max = iC^^ K, at an observing frequency of 
1 GHz, so that 7 « 200, then, from Eq. ([T5]), the number 
of terms contributing to the sum is A^max = 2. 
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4. Stationary solutions 

4.1. Intra-day variable sources 

For comparison with observations of intra-day variable 
sources, it is convenient to formulate the expression for 
the specific intensity given in Eq. (IA.15|) in terms of quan- 
tities accessible to observation. Expressing the result in 
terms of the observed (at z — 0) brightness temperature, 
we find, in the case of weak absorption, and at low fre- 
quency (i/ ^ I/g) 



(1 + .)' 



1/5 



.1/5 



2/15 -1/3 ,r,^\ 
^^maxM'^GHz K (27) 



(jKirk fc Tsand l2006f ) where V = lOPio is the Doppler 
boosting factor, z is the redshift of the host galaxy, Ts is the 
optical depth of the source at the observing frequency v = 
VGUz GHz, and the characteristic synchrotron frequency 
of the electrons is i^s — ^'inax,i4 x 10^^ Hz. 

According to Eq. ([27| . brightness temperatures 
of Tb « 10^'^ K, such as observed in the sources 
PKS 1519 - 273 and P KS 0405 -385 (iMacauart et al 



20001: iRickett et al.ll2002t ) can be understood within a sim- 
plified homogeneous synchrotron model in which ^ ^ 1, 
implying a relatively modest inverse Compton luminosity, 
i.e., no catastrophe. Even the extremely compact source 
J 1819 -h3845, which has Tb ^ 2 x lO^^^K can be ac- 
commodated in a catastrophe-free model provided the 
Doppler factor is greater than about 15. In each case, a 
hard spectrum is predicted, extending to fmax,i4 x 10^^ 
Hz. Although the dependence of the brightness temper- 
ature on this parameter is quite weak, si multaneous ob 



serva tions in the radio to IR and optical (jOstorero et al 



200i) have the potential to rule out this explanation on a 
source by source basis. 

4.2. Resolved sources 



In a seminal paper. iReadheadl ([199J) discussed the distri- 
bution in brightness temperature of a sample of powerful 
sources whose angular size could either be measured di- 
rectly, or constrained by interplanetary scintillation. In 
discussing these objects several simplifications must be 
made, even within the context of a homogeneous syn- 
chrotron model. 

Firstly, in the two low fr equency samples (81.5 MHz 
and 430 MHz) considered bv IReadheadl (|l994f) . the emis- 
sion is thought to be almost isotropic. Doppler boosting is 
then unimportant and can be neglected. Secondly, these 
sources are not very compact; their extension on the sky 
is typically between 0.1 and 1 arcsec. Therefore, for our 
discussion we fix the linear extent R of the source to 
1 kpc, corresponding to an angular size of approximately 
0.2 arcsec at redshift z — 1. This leaves three parameters 
needed to specify the source model: the magnetic field 
strength B, the electron density rio and the Lorentz fac- 
tor 7 of the electrons. In order to clarify the physics of 
a source, we transform from the parameter set {B,ne,j) 



to the characteristic Lorentz factors 7cq and 7cat defined 
in Eqs. ((23|) and ([24l) . Our basic parameter set is there- 
fore (7oq, 7cat, 7)- Finally, in order to display on a two- 
dimensional figure source properties such as brightness 
temperature and spectral slope at a particular frequency, 
we consider a slice through this three dimensional param- 
eter space, selecting parameters such that the particle and 
magnetic energy densities are in equipartition: 7 = 7oq. 

The properties of source models on this slice are shown 
in the 7eq-7cat plane in Fig. [2l This plane can immediately 
be divided into regions of strong and weak absorption, as 
defined in Eq. ^ . The boundary, drawn as a thick dashed 
line, represents the locus of the points at which 7c — 7cq. 
Weakly absorbed sources lie towards higher 7cq and 7cat 
(i.e., the upper-right side) and strongly absorbed sources 
towards lower 7eq and 7cat (i-e., the lower-left side). We 
also show (in blue) contours of the magnetic field strength. 

The remaining source properties depend upon the 
choice of observing frequency. In Fig. [2] we take this to 
be 81.5 MHz , corresponding t o the low frequency sample 
discussed bv IReadheadl (|l994l ). In order to determine the 
spectral slope of a given source, we plot as a yellow line 
the locus of points where the observing frequency coin- 
cides with the frequency at which the optical depth to 
absorption is unity, Uahs- Sources that lie above this line 
(on the side of larger 7cat) are optically thin at the cho- 
sen observing frequency. In addition, the red line in Fig. [2] 
gives the locus of points where the observing frequency 
equals the characteristic frequency of synchrotron radia- 
tion Vc- By definition, the intersection point of these lines 
lies also on the boundary between weak and strong ab- 
sorption (the thick dashed line) . The observing frequency 
lies below I'c on the lower-right side of the red line. The 
colour shading gives the intrinsic brightness temperature 
at the chosen observing frequency. 

The two lines (yellow and red) divide the 7oq-7cat-pla'ne 
in Fig. [2] into four regions with differing spectral proper- 
ties: in region A, sources have a Rayleigh- Jeans spectrum 
Ii, cx in region B, the spectrum is that of low fre- 
quency, optically thin synchrotron ra diation I,^ cx v^^^, in 
region C, it is close to (x v (see Slvsh 19921 ) and in 
region D it falls off exponentially Ii, oc exp(— i//i/c). 

Consequently, flat spectrum sources reside in region B, 
preferentially close to the yellow line and in region C, pref- 
erentially close to the red line. 

Sources that are in equipartition and lie below the 
threshold of the Compton catastrophe are to be found in 
the upper left half of Fig.|2l above the dotted line on which 
7oq = 7cat- The maximum brightness temperature accessi- 
ble to these sources occurs close to 7cq = 7cat = 10^, and 
is approx imately 10^^'^ K, in rough a greement with the 
results of iKellermann fc Pauhnv-TothI ([1969), who, how- 
ever, did not assume their sources to be in equipartition. 
The brightest sources are weakly absorbed, (they lie to 
the right of the thick dashed line) and have a magnetic 
field strength of a few milliCauss. Their optical depth to 
synchrotron self-absorption lies close to unity at the ob- 
servation frequency (they lie close to the yellow line). 
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Fig. 2. The brightness temperature as a function of 7cq 
and 7cat assuming equipartition between the magnetic and 
particle energy densities and a source size 1 kpc. Black 
contour lines indicate logio(T/Kelvin) ~ 9, 10, 11, 12, 
13 and 14. The red line is the locus of points at which 
the characteristic synchrotron frequency of the emitting 
particles is 81.5 MHz, the yellow line shows where the 
source has an optical depth of unity at this frequency. 
The dashed line divides regions of strong absorption (to 
the left) from those of weak absorption (to the right). The 
diagonal 7cq = 7cat is shown as a dotted line. Contour lines 
of the magnetic field strength are shown in blue, ranging 
from log]^Q(_B/Gauss) = —4 to (in the bottom right-hand 
corner). 



Sineal fc Gopal-Krishnal (Il985l) first discussed the ef- 



fects of the additional assumption of equipartition on 
bright sources and used it to esti mate Dopp l er fac tors 



for rapidly variable sources. Later, iReadheadI (|1994D in- 
troduced the concept of an "equipartition brightness tem- 
perature" to explain the observation that the temperature 
distribution of resolved sources appears to peak signifi- 
cantly below 10^^ K. However, the crucial additional as- 
sumptions in his treatment is that the source flux is mea- 
sured at the "synchrotron peak", and that the electron 
distribution is a power-law in energy. This implies that 
the opacity at a given frequency (e.g., at the synchrotron 
peak) is dominated by those electrons with a correspond- 
ing characteristic frequency. In our model, in which the 
electron distribution is approximated as monoenergetic, 
these assumptions are roughly equivalent to demanding 



that the source lies on the red line in Fig. [2] if it is weakly 
absorbed (i.e., on the boundary of regions B and D), and 
on the yellow line if it is strongly absorbed (i.e., on the 
boundary of regions C and D). This leads to a maximum 
bri ghtness te r npera ture of a few times 10^° K , as found 
bv 'R eadheadI (1994). Furthermore, as noted bv lReadheadI 
(1994), such sources lie far from the threshold tempera- 
ture, achieved along the dotted line in Fig. [21 

Replacing the assumption that the source flux is 
measured at the synchrotron peak, by the requirement 
that its spectrum be flat, i.e., that it lie in region B 
of Fig. [21 one sees that a wide range of brightness 
temperatures is available for sources in equipartition, 
extending up to the thres h old t emperature found by 
Kellermann fc Paulinv-Toth ( 1969( ). Thus, the observed 



temperature distribution is not explained by the assump- 
tion of equipartition. 

5. Time dependence and acceleration 

In order to explain the occurren ce of brightness temper- 
atures above 10^^ K, ISlyshI (|l992l ) formulated a model in- 
volving a monoenergetic electron distribution in a strongly 
absorbed source, in the sense that 7 < 7c, where 7c is de- 
fined in Eqs. ^ and ([3]). He considered two scenarios, (i) 
a time-dependent one in which electrons were injected at 
arbitrarily high Lorentz factors and allowed to cool and 
(ii) one in which a strong continuous re-acceleration of 
the electrons led to a high brightness temperature equi- 
librium. 

In each case, the assumption that the source is strongly 
absorbed leads to extreme values of the parameters. For 
example, in the first scenario in which high en ergy parti- 
cles are injected into the source, ISlyshI (jl992i ) finds that 
a brightness temperature of Tb > 5 x 10^^ K can be sus- 
tained over 1 day at an observing frequency of 1 GHz. 
This is clearly in conflict with our analysis. The elec- 
tron Lorentz factor required to achieve this temperature is 
7 > 10^. However, the condition that the source is strongly 
absorbed, which is used in this model to estimate the cool- 
ing rate, combined with the condition i/g « 1 GHz required 
for a flat spectrum, leads to an extremely large Thomson 
optical depth, t « 130, as well as an implausibly low mag- 
netic field B Ri 2 X 10~^^ G. The parameter ^ that deter- 
mines the inverse Compton luminosity is approximately 
10^^, which implies an extremely large compactness of the 
inverse Compton radiation from the source. The result- 
ing copious pair production invalidates the analysis and, 
ultimately, reduces the brightness temperature achievable 
in the radio range. The same criticism app lies also to the 
second scenario described by ISlvshI (Il992i ) in which ac- 
celeration balances inverse Compton losses to provide a 
brightness temperature of lO^'^K at 1 GHz. 

In the absence of Klein-Nishina effects on the scatter- 
ing cross section, we find the time dependence of the parti- 
cle and photon energies can be described by the three ordi- 
nary differential equations ^U^ . (|^T|) and ((221) ■ Inspection 
of these shows that if the threshold temperature is ex- 
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ceeded (7 > 7cat)i the inverse Compton luminosity grows 
in a timescale of roughly the light-crossing time of the 
source. Thus, the threshold can only be substantially ex- 
ceeded if the acceleration process in Eq. operates on a 
shorter timescale. However, these equations employ a spa- 
tial average over the emission region. Although a rapid 
acceleration rate might be achieved locally in small re- 
gions of the source, once an average is taken, no timescale 
in the system can be shorter than the light-crossing time 
of the region over which the accelerated particles are dis- 
tributed. In this case, the threshold temperature cannot 
be significantly exceeded. 

At first sight, Klein-Nishina effects offer a possible es- 
cape from this conclusion. If even the first order scattering 
is suppressed, which requires extremely large Lorentz fac- 
tors for the electrons (7 > 10^° is needed for Klein-Nishina 
effects when scattering 10 GHz photons), the strong reduc- 
tion in the rate of cooling by inverse Compton scattering 
suggests that higher brightness temperatures Tb might be 
possible. 

This is, however, not the case, because the rate of pro- 
duction of electron-positron pairs by photon-photon in- 
teractions becomes important. The strength of this effect, 
which is not included in our model equations, can be mea- 
sured in terms of the "compactness parameter" i (see, for 
example Mastichiadis fc Kirk 1995f) . defined as 

i = II^Elh^ (28) 

where Un„,^^ is defined in Eq. A^max in Eq. (fT5|) . 

and vn^^^ is taken to be (47^/3)^"''''i'o- When i > \, 
one expects the pair-production rate to be roughly equal 
to the light-crossing time of the source. This leads to a 
sharp rise in the Thomson optical depth, invalidating the 
assumption of scatter-free escape of synchrotron photons 
that is implicit in our model. The associated confinement 
of these photons reduces the brightness temperature. 

We illustrate this in Fig. [31 where we compare two 
models with the same linear size R (and observing fre- 
quency), but different electron densities and different 
values of B, chosen as follows: For any given set of pa- 
rameters, R, B and ric, and observing frequency fobs, 
the optical depth to synchrotron absorption Tg, as de- 
fined in Eq. (jA.12p . has a single maximum as a func- 
tion of 7, located close to the point where fobs equals 
the characteristic synchrotron frequency. If the source is 
optically thick to absorption at this point, then 7 < 7c, 
as described in Sect.[2l and the brightness temperature is 
roughly Z^mc^ / Ak-Q- If, on the other hand, the source is 
optically thin at this point, then 7 > 7c, but the brightness 
temperature, given approximately by Tg x 37mc^/4fcB, de- 
creases to higher 7, as can be seen from Eq. (|A.12[) . Thus, 
assuming inverse Compton scattering does not intervene, 
the maximum brightness temperature is observed at a fre- 
quency such that Tg « 1, when 7 = 7c, which implies 
a; « 1. These conditions are imposed on the parameters of 
the models presented in Fig. [31 In addition to the source 
size, chosen to be i? = 0.01 pc and the observing frequency. 
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Fig. 3. Upper panel: The brightness temperature Tb 
(blue), and the compactness £ (red) as functions of time, 
for two stationary, local sources (I? = 1, z = 0) with lin- 
ear size R = 0.01 pc, observed at 1 GHz. The Thomson 
optical depth is tt = 0.01 (dashed lines) and tt = 1 (solid 
lines) and the remaining parameters are chosen such that 
the optical depth to synchrotron self- absorption Tg ~ 1 at 
7 = 7c (see Eqs. © and A horizontal line is drawn 
to indicate I — 1. Lower panel: The electron Lorentz fac- 
tor (purple dashed) and the optical depth to synchrotron 
self-absorption Tg (green dashed) for the case tt = 0.01. 
A horizontal line indicates Tg = 1. 



set to 1 GHz, this leaves one free parameter, which we 
choose to be the optical depth to Thomson scattering tt . 

The upper panel in Fig. [31 shows the time-dependence 
of the brightness temperature found by solving Eqs. (|12p 
and (fTT]) numerically for sources with tt — 0.01 (dashed 
blue line) and tt = 1 (solid blue line), without al- 
lowance for Doppler boosting {T> = 1). These sources 
have 7c = lO"^'^ and 7c = 10* '^, respectively and, in the 
absence of inverse Compton cooling, they could poten- 
tially achieve brightness temperatures of Tb « lO^'^^'^K 
and Tb w lO^'^'^K. In order to do so, rapid accelera- 
tion is required, since for these source parameters, inverse 
Compton cooling leads to a time-asymptotic value of the 
Lorentz factor that is somewhat lower than 7c for slow 
acceleration. The exact value of the asymptotic solution 
depends on the strength of the acceleration. For accel- 
eration on the light-crossing timescale, it corresponds to 
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a ~ 7/7tr [see Eq. (f^ ]. In Fig.[3]we choose a = 1.57/7tr, 
which leads to an overshoot that shghtly exceeds 7c. 

For tt = 0.01, the compactness, shown as a function 
of time by the red dashed line, remains well below unity, 
so that the effects of pair production can be neglected. 
However, this is not the case for tt = 1. Here, the com- 
pactness (solid red line) rises rapidly, reaching unity at 
i « 0.25, where Tb ~ 3.5 x 10^^ K, well below its potential 
maximum. Thus, the attempt to gain higher brightness 
temperature by increasing tt, and, hence, 7c, leads to a 
breakdown in the model assumptions due to pair produc- 
tion. 

The lower panel of Fig. [3] shows the electron Lorentz 
factor and the optical depth to synchrotron self-absorption 
Ts as functions of time for the case tt = 0.01. The Lorentz 
factor overshoots both its time-asymptotic value and 7c. 
Correspondingly, the optical depth, (shown as the green 
dashed line) which initially rises with 7, reaching unity 
at 7 = 7c goes through a maximum very shortly after- 
wards. However, the overshoot is not sufficient to push Tg 
back below unity, and the maximum brightness tempera- 
ture, which coincides with the maximum Lorentz factor, 
remains at Tb = 5 x 10^^ K, somewhat below the value of 
Tb « lO^'^-^ K, estimated for large optical depth. 

6. Conclusions 

The well-known upper limit on the brightness tempera- 
ture of a synchrotron source Tb <^ 10^^ K imposed by the 
inverse Compton catastrophe, has been reassessed, assum- 
ing a monoenergetic electron distribution. 

In weakly absorbed sources (see Eq. ([2])), this distri- 
bution mimics the situation in which the conventional 
power-law is truncated to lower energies at a Lorentz fac- 
tor 7niin. Using the standard theory of synchrotron emis- 
sion and self-absorption, we find that, for such sources, 
the brightness temperature at a frequency of a few GHz 
can reach approximately 10^^ K, the precise limit being 
given in Eq. ([27|) . Physically, this increased limit reflects 
the absence of cool electrons in monoenergetic distribu- 
tions and in those that are truncated or hard below a cer- 
tain Lorentz factor. As a consequence, intra-day variable 
sources can in principle be understood without recourse 
to ot her mechan isms such as unusually la rge Doppler fac- 
tors coherent emission (e. g.. iBegelman et al 



2005f ) or proton synchrotron radiation ( KardashevI 2000[ ) 



The possibility of exceeding the new limit in a time- 
dependent solution by balancing losses against a strong 
acceleration term has been investigated using a set of spa- 
tially averaged equations. Provided the acceleration pro- 
cess remains causal i.e., the acceleration time averaged 
over the source remains longer than the light-crossing 
time, we find a modest overshoot is possible, but the 
maximum temperature is still restricted by Eq. (I27|) . In 
strongly ab sorbed sources, such as those considered by 
Slvshl (|l992h . high brightness temperatures cannot be at- 



tained in a self-consistent model of the kind we discuss. 
The underlying reason is that extremely compact sources 



would be required, in which copious pair-production must 
be taken into account. 

We have examined in detail the parameter space avail- 
able to homogeneous synchrotron sources of fixed size. In 
the case of flat spectrum sources, we find that the impo- 
sition of the condition of equipartition between the par- 
ticle and magnetic field energy densities does not result 
in a lower limit on the brightness temperature than that 
given by the in verse C ompton catastrophe. Suggestions 
to the contrary (jReadh ead 1994) are based on the more 
restrictive twin assumptions that the power-law electron 
distribution is not truncated within the relevant range, 
and that the temperature is measured at the point where 
the optical depth of the source is approximately unity. 
Consequently, the observed temperature distribution does 
not support the equipartition hypothesis. We also find 
that flat spectrum sources close to equipartition can ap- 
proach the threshold temperature of the inverse Compton 
catastrophe, in contrast with the finding based on the 
more restrictive assumptions in iReadheadI (|l994l ). 

A corollary of the theory presented here is that for 
very bright sources, which are necessarily weakly ab- 
sorbed, a high degree of intrin sic circular polarisation is 
predicted ( Kirk fc Tsand 2006f) . In addition, the theory 
can be tested by comparison of the predicted synchrotron 
spectrum with simultaneous observations of high bright- 
ness temperature sources in the radio to infra-red range 
and comparison of the predicted inverse Compton emis- 
sion with measurements at MeV to GeV energies. 



Appendix A: Synchrotron formulae 

We consider a region of homogeneous magnetic field S, lin- 
ear dimension R, (and volume R^) containing monoener- 
getic electrons/positron of number density Ue and Lorentz 
factor 7. The volume emissivity for synchrotron radiation, 
summed over polarisations, is: 



X 



V3 



n^af fiVL^ sin 9 F {x) 



47r 

ar^L sin 6*72 



4Tr 

VQl^ K = 3f^Lsin6'/(47r)] 



F{x) = X / dii^5/3(0 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



where ttf = jfic is the fine-structure constant, JIl = 
eB /mc the Larmor frequency and 9 the angle between the 
magnetic field and the direction of the emitted radiation. 
For small and large x the limiting forms are 



F{x) 



in 



(- 

V3r(l/3) V2 



1/3 



for a; <C 1 



F{x) a/~2~*^ ^ for a; ^ 00 



(A.5) 
(A.6) 
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The absorption coefBcient for unpolarised radiation is 

1 UeaT Kq/^{x) 



2\/3Q!f6sin^ 7^ 



To find the energy density Us in synchrotron photons 
in a given source, I^, must be integrated over angles and 
(A. 7) over frequency. The result depends on the geometry and 
optical depth as well as the position within the source. 



where h = n^^/m(? is the magnetic field in units of the However, an average value can be estimated by introduc- 
critical field = ^A\A x G and ax is the Thomson ^ geometry dependent factor C « 1: 



cross-section. The limiting forms are: 
K„^{x) « ^!^!£g/^ fora;«l 



3-5/3 

TT 



for 



(A.8) 
(A.9) 



c Jo 



(A.19) 



and denoting by {ly) the specific intensity evaluated at 
e = tt/2. Then 



Because a^, is a monotonically decreasing function of 
X, we can define a unique x^ib, 7) where the optical depth 
Ts = Rav for synchrotron absorption along a path of 
length R is unity: 



(A.20) 



with 



Rai, (xa) = 1 



(A.IO) 



roo 

C/(7) = 7' / da;5(7,x){l-exp[-r,(7,x)]} (A.21) 
Jo 



This integral is dominated by the region x ^ x^, in. the 
If Xa <C 1, we have weak absorption and for Xa » 1 strong .^yg^k absorption regime: 
absorption. The transition between the two regimes occurs 
near Lorentz factor 70, defined as 



C/(7) 



7c = 
so that 



2V3af6 sin( 



1/5 



f 

Jo 



dx Sts 



(A.11) 



dxF{x) 



Ts = 1 ^ if 5/3 (a;) 

_ a/Stttoc^ if 5/3(3^) 
87re^fc7^ 

where 7 = 7/70- In the case of weak absorption, 

Xa « 2^/5 [r(5/3)]'/' /f for 7 » 1 
whereas in the strong absorption regime 
Xa ~ — 51n7 for 7 <C 1 



for 7 > 1 



(A.22) 



and by the region around x = Xa in the strong absorption 
(A.12) regime: 

f/(7) « / dxS' 
Jo 

(A.13) « 12.57^ (In 7)^ for7<l (A.23) 

which suggests the simple approximation 
12.57^ 



(A.14) 



f/(7) 



0.183+ (In 7)^ 



(A.24) 



7.7575 



Taking account of synchrotron emission and absorp)- 
tion and ignoring the role of polarisation, the specific in- 
tensity in a direction which cuts the source on a path of ^^^^.^ *he constant 0.183 was chosen such that the ap 



length R is 

= S^[l- exp {-Ts)] 
where the source function Si, = jv/ctv is 



proximation passes through the point U{1) = 0.945 found 
by numerical integration. 
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